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Abstract 

Aircraft  pitch  response  is  a  crucial  element  of  piloted  vehicle  flying  qualities.  The 
short  term  pitch  response  has  created  controversy  over  the  form  and  substance  of  any 
requirements.  Currently  there  are  six  different  methods  for  evaluation  in  MIL-STD-1797A. 
There  are  many  other  methods  which  have  been  proposed.  The  biggest  problem  is  that 
many  of  these  methods  often  give  conflicting  results.  The  overall  goal  of  the  present  effort 
is  to  compare  and  contrast  the  Time  Domain  criterion  and  the  Ralph  Smith  criterion.  By 
examining  these  methods  on  common  grounds,  areas  of  agreement  and  discrepancies  can  be 
found.  Parametric  studies  are  performed  and  trends  identified. 


COMPARISON  OF  THE  RALPH  SMITH  AND  THE  TIME 
DOMAIN  FLYING  QUALITIES  CRITERIA 


1.  Introduction 

There  are  many  elements  in  aircraft  flight  mechanics  which  contribute  to  the  overall  flying 
qualities  of  an  aircraft.  These  include  the  pitch  response,  roll  response,  and  yaw  response. 
One  of  the  most  important  is  the  short-term  pitch  response  of  an  aircraft.  MIL-STD- 
1797A,  Flying  Qualities  of  Piloted  Aircraft  [1],  offers  six  different  methods  for  evaluating 
short-term  pitch  response.  Each  method  has  strengths  and  weaknesses  depending  on  aircraft 
classification  and  flight  phase.  All  six  methods  have  been  maintained  because  the  short¬ 
term  pitch  response  characteristics  are  regarded  as  important[l].  MIL-STD-1797A  provides 
some  guidance  for  determining  the  appropriate  method  to  apply.  Still,  one  must  decide 
upon  which  of  the  six  methods  to  include  in  a  speciflcation.  The  simple  answer  of  including 
all  six  will  lead  to  conflicting  results  [1].  There  are  also  recognized  methods  that  are  not  in 
MIL-STD-1797A,  but  are  used  by  people  in  the  industry,  such  as  Numerator  Time  Constant, 
Bandwidth  and  Phase  Sensitivity,  and  the  Ralph  Smith  criterion  [2]  that  create  similar  types 
of  conflict. 

Two  methods  are  analyzed  in  this  research  -  the  Time  Domain  criterion,  from  MIL- 
STD-1797A  [1],  and  the  Ralph  Smith  criterion,  from  AFFDL-TR-78-154  [3].  By  comparing 
these  criteria,  it  can  be  seen  where  regions  of  conflict  and  agreement  are  located.  This 
will  aid  in  making  a  decision  to  either  keep  the  criteria  the  way  they  are,  modify  them,  or 
possibly  even  combine  them  into  one  new  criterion. 

1.1  Overview 

The  study  of  flying  qualities  is  the  discipline  in  aeronautical  engineering  that  is  concerned 
with  basic  aircraft  stability  and  controllability.  ‘Flying  qualities,’  ‘stability  and  control,’ 
and  ‘handling  qualities’  are  three  terms  which  are  generally  considered  synonymous  [4].  To 
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Figure  1.1:  Flying  Qualities  Breakdown 


prevent  confusion,  the  following  definitions  will  be  used.  Both  the  US  Air  Force  Test  Pilot 
School  and  US  Naval  Test  Pilot  School  agree  that  ^‘'flying  qualities  are  those  stability  and 
control  characteristics  which  influence  the  ease  of  safely  flying  an  aircraft  during  steady 
and  maneuvering  flight  in  the  execution  of  the  total  mission”  [5] .  Edkin  defines  stability  as 
“...the  tendency  or  lack  of  it,  of  an  airplane  to  fly  with  wings  level”  and  control  as  “...steering 
an  airplane  on  an  arbitrary  flight  path”  [6].  Cooper  and  Harper  define  handling  qualities 
as  “...those  qualities  or  characteristics  of  an  aircraft  that  govern  the  ease  and  precision 
with  which  a  pilot  is  able  to  perform  the  tasks  required  in  support  of  an  aircraft  role”  [7]. 
Fig.  (1.1)  shows  how  flying  qualities,  stability,  control,  and  handling  qualities  are  related  to 
one  another.  The  figure  shows  that  stability,  control  and  handling  qualities  are  all  subgroups 
of  flying  qualities. 

Stability  and  control  analysis  deals  with  the  interaction  of  the  control  surfaces  with  the 
external  forces  and  moments  on  the  aircraft  [8].  Generally,  stability  and  control  analysis 
primarily  deals  with  systems  that  are  still  in  the  design  phase,  as  shown  in  Fig.  (1.2).  Note 
that  there  is  no  pilot  in  this  system.  A  stick  force,  Fg,  is  applied,  which  goes  through  the 
control  system,  to  produce  an  elevator  deflection,  Sg.  This  elevator  deflection  is  fed  into  the 
aircraft  dynamics  to  produce  the  desired  pitch  response. 


_ Fs 

Control 

Aircraft 

0 

System 

Dynamics 

Figure  1.2:  Open  Loop  System 
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On  the  other  hand,  handling  qualities  assessment  deals  with  the  pilot  and  aircraft  performing 
as  a  closed-loop  system,  as  shown  in  Fig. (1.3).  The  pilot  wants  a  desired  pitch  angle,  0,  so 
he  puts  in  a  stick  force.  The  input  goes  through  the  control  system  and  aircraft  dynamics. 
The  resulting  pitch  angle  is  feed  to  the  stability  augmentation  system,  SAS,  and  back  to 
the  pilot,  where  he  determines  if  more  input  is  needed.  The  main  difference,  then,  is  that 
stability  and  control  require  analysis  without  a  pilot,  while  handling  qualities  is  considered 
an  analysis  with  a  pilot  in  the  loop.  In  this  research,  the  term  open-loop  is  used  to  signify 
that  there  is  no  pilot  in  the  analysis  of  the  aircraft  transfer  function,  while  closed-loop  means 
a  pilot  model  is  in  the  analysis. 


Figure  1.3:  Closed  Loop  System 


The  six  methods  offered  by  MIL-STD-1797A  are  used  to  predict  handling  qualities  while 
the  aircraft  is  still  in  the  development  phase.  Since  aircraft  perform  a  wide  variety  of 
maneuvers  and  vary  in  size,  some  type  of  grouping  is  necessary  before  analysis  can  begin. 
The  class  designations,  as  described  in  MIL-STD-1797A,  are  used  to  help  determine  the 
requirements  according  to  broad  categories  of  intended  use.  The  intended  use  of  an  aircraft 
must  be  known  before  required  configurations,  loadings  and  operational  flight  envelopes 
can  be  defined.  Four  classes  of  aircraft  are  defined  by  MIL-STD-1797A  and  are  outlined  in 
Table  1.1. 

This  research  examines  a  range  of  values  representative  of  two  different  aircraft.  The 
first  is  Calspan  Corporation’s  variable  stability  Learjet  24  shown  in  Fig.  (1.4).  The  second  is 
the  Variable  Stability  In-flight  Simulator  (VISTA)  which  uses  the  F-16D  as  its  host  aircraft, 
pictured  in  Fig.(1.5).  The  Learjet  can  be  considered  a  Class  II  aircraft,  while  the  F-16 
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Table  1.1:  Aircraft  Classification 


Class 

Description 

Example 

I 

small  light  aircraft 

T-41,  OV-10 

II 

medium  weight  aircraft  with 
low-to-medium  maneuverability 

C-21,  C-130 

III 

large,  heavy  aircraft  with 
low-to-medium  maneuverability 

KC-10,  B-2 

IV 

highly  maneuverable  aircraft 

F-16,  F-117 

VISTA  can  be  considered  a  Class  IV  aircraft.  Using  variable  stability  aircraft  makes  it 
possible  to  flight  test  the  analytic  results  of  this  present  study  in  future  research.  The 
cockpit  environment  can  be  changed  to  match  that  of  another  aircraft.  Since  an  aircraft  is 
being  flown,  it  provides  a  degree  of  realism  that  cannot  be  duplicated  in  a  ground-based 
simulation.  As  the  pilot  moves  the  controls,  he  can  experience  the  true  flight  motions, 
accelerations,  and  handling  qualities  of  the  simulated  aircraft  [9].  This  realism  gives  the 
pilot  a  higher  level  of  confidence  when  determining  a  handling  qualities  level. 

In  this  research,  a  mapping  will  be  provided  to  show  regions  of  agreement  and  conflict 
for  two  of  the  different  methods  used  to  predict  handling  qualities.  Since  every  point  in  the 
region  represents  a  different  aircraft  transfer  function,  a  variable  stability  aircraft,  such  as 
the  VISTA  .F-16,  can  be  configured  to  represent  one  of  these  transfer  functions.  A  point 
in  a  region  of  conflict  could  be  programmed  into  the  flight  control  system  so  a  test  pilot 
could  fly  the  simulated  aircraft  and  determine  the  real  handling  qualities  level.  If  enough 
cases  are  flight  tested  and  the  results  show  a  specific  trend,  a  decision  could  be  made  on 
whether  one  of  the  methods  used  to  predict  handling  qualities  needs  to  be  modified  or  even 
eliminated. 

Experience  with  aircraft  operations  indicate  that  certain  flight  phases  require  more  strin¬ 
gent  values  of  flying  qualities  parameters  [1].  MIL-STD-1797A  defines  three  categories  of 
flight  phases,  outlined  in  Table  1.2.  This  research  only  examines  the  Category  C  flight 
phase,  approach  and  landing.  However,  the  method  derived  in  Chapter  III  of  this  present 
eflPort  can  handle  any  flight  phase  category.  With  aircraft  classification  and  flight  phase 
known,  handling  qualities  levels  can  be  addressed. 
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Figure  1.4:  Lear  Jet 


Figure  1.5:  F-16  VISTA 


Table  1.2:  Flight  Phase  Categories 


Category 

Description 

Example 

A 

nonterminal  flight  phases  that  require 
rapid  maneuvering,  precision  tracking, 
or  precise  flight-path  control 

air-to-air 

combat 

B 

nonterminal  flight  phases  that  require 
gradual  maneuvering  without 
precision  tracking 

climb 

C 

terminal  flight  phases  that  require 
gradual  maneuvering  with 
accurate  flight-path  control 

landing 

When  determining  the  handling  qualities  of  an  aircraft,  a  pilot  must  answer  a  series  of 
questions.  These  questions  lead  to  a  pilot  opinion  rating  on  the  Cooper-Harper  scale  [7]. 
The  scale  was  developed  over  the  years  as  a  means  of  putting  a  short-hand  symbol  to  the 
comments  the  pilots  made  about  a  aircraft  performance  [7].  By  answering  a  series  of  yes-no 
questions,  the  pilot  can  relate  the  controllability,  the  workload  requirements  and  the  amount 
of  improvement  needed. 

MIL-STD-1797A  defines  three  levels  of  handling  qualities,  outlined  in  Table  1.3  [1]. 
These  levels  are  based  on  the  the  Cooper-Harper  scale,  shown  in  Fig.  (1.6)  [1].  A  correlation 
between  the  Cooper-Harper  scale  and  the  handling  qualities  levels  defined  by  MIL-STD- 
1797A.  A  Cooper-Harper  rating  of  1  to  3  defines  the  Level  1,  4  to  6  defines  Level  2,  and  7 
through  9  are  Level  3  [7]. 

1.2  Previous  Work 

Work  done  on  the  topic  of  handling  qualities  comparison  is  not  limited  to  this  research. 
Research  is  done  on  a  continuing  basis  at  the  Flight  Dynamics  Directorate  of  Wright  Labo¬ 
ratory.  The  Handling  Qualities  group  is  currently  working  on  comparison  mappings  for  the 
different  flight  phases.  This  work  is  on  going  in  an  attempt  to  determine  which  criteria  to 
include  in  a  revision  of  MIL-STD-1797A. 

Another  source  of  work  done  was  performed  by  Kish  [10]  for  his  Master’s  Thesis  at 
the  Air  Force  Institute  of  Technology.  His  research  examined  the  Neal-Smith  criterion  and 
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Pilot  compensation  not  a  factor  for 
desired  peiformance 

1 

Pilot  compensation  not  a  factor  for 
desired  performance 

2 

Minimal  pilot  compensation  required  for 
desired  penormance 

3  ^ 

Desired  performance  requires  moderate 
piFot  cornpensation 

4 

Adequate, performance  requires 
con^derame  pilot  compensation 

5 

Adequate  performance  requires 
extensive  pilot  compensation 

6 

Adequate  performance  not  attainable  with 
maximum  tolerable  pilot  compensation. 
Controllability  not  in  question 

N 

7 

Considerable  pilot  compensation  is  required 
for  control 

8 

Intense  pilot  .compensation  is  required  to 
retain  control 

% 

Control  will  be  lost  during  some  portion  of 
required  operation 


10 


S  DESIGNATION  OF  FLIGHT  PHASE 
)mONS 


Scale 


Table  1.3:  Handling  Qualities  Levels 


Level 

Meaning 

Description 

1 

satisfactory 

Flying  qualities  clearly  adequate  for  the  mission 
flight  phase.  Desired  performance  is  achievable 
with  no  more  than  minimal  pilot  compensation. 

2 

acceptable 

Flying  qualities  adequate  to  accomplish  the  mission 
flight  phase,  but  some  increase  in  pilot  workload 
and/or  degradation  in  mission  effectiveness  exists. 

3 

controllable 

Flying  qualities  are  such  that  the  aircraft  can  be 
controlled  in  the  context  of  the  mission  flight 
phase,  even  though  pilot  workload  is  excessive 
and/or  mission  effectiveness  is  inadequate. 

the  Usp,Csp:T0^,Te  criterion.  The  Neal-Smith  criterion,  as  well  as  the  less  restrictive  Pilot- 
in-the-Loop  criterion,  uses  an  optimal  pilot  model  to  predict  the  handling  qualities  of  an 
aircraft.  Kish  compares  both  of  these  closed  loop  methods  with  the  open  loop  ujgp,  Csp, 
criterion,  in  order  to  compare  the  results,  in  the  landing  phase  of  flight. 

Work  on  yet  another  topic,  similar  in  nature,  is  also  being  flight  tested  at  the  US  Air 
Force  Test  Pilot  School,  by  Kivioja  [11].  In  that  research,  the  control  anticipation  parameter, 
CAP,  is  compared  to  the  Bandwidth  the  landing  phase.  In  the  test  phase  of  the  research, 
pilots  will  evaluate  the  landing  task  using  the  Cooper-Harper  rating  scale  [11].  These  results 
will  then  be  compared  to  the  predicted  ratings  received  by  evaluation  of  the  criteria. 

1.3  Research  Objectives 

The  overall  goal  of  the  present  effort  is  to  compare  and  contrast  the  Time  Domain 
criterion  and  the  Ralph  Smith  criterion.  In  order  to  accomplish  this  overall  goal,  some 
specific  objectives  must  be  met. 

1.  Develop  a  computer  code  enabling  the  two  criteria  to  be  examined  at  each  transfer 
function  in  a  grid  of  the  undamped  natural  frequency,  Ugp  versus  the  damping  ratio,  (sp- 

2.  Develop  a  mapping  system  to  identify  the  areas  of  conflict  and  agreement  between  the 
criteria. 

3.  Perform  a  parametric  study  to  identify  trends  for  areas  of  agreement. 
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1.4  General  Approach 


An  aircraft  transfer  function  can  be  described  by  five  parameters:  uJsp,T82tCsp->K9  and 
Tg.  By  holding  Tg^,  Kg  and  rg  constant  while  varying  cogp  and  Csp  a  two  dimensional  region 
of  aircraft  transfer  functions  can  be  created.  Once  each  transfer  function  is  determined, 
the  handling  qualities  level  according  to  the  Ralph  Smith  criteria  is  then  determined.  This 
creates  a  map  in  a  grid  of  ojgp  vs.  (sp  of  the  Ralph  Smith  criteria.  The  same  procedure  is 
done  for  the  Time  Domain  criteria.  The  two  criteria  are  then  compared  to  determine  the 
areas  of  conflict  and  agreement. 

1.5  Overview 

This  research  is  separated  into  four  chapters.  Chapter  If  contains  the  background 
information  necessary  for  understanding  the  short-period  pitch  response  of  an  aircraft.  Also 
included  is  the  background  material  of  the  criteria  involved  in  the  mappings.  Chapter  III 
describes  the  results  of  completing  Objective  1  and  Objective  2.  A  sample  mapping  is 
provided  to  illustrate  the  algorithm.  The  results  of  the  parametric  study,  Objective  3, 
is  presented  in  Chapter  IV.  Chapter  V  gives  the  conclusions  which  were  drawn  from  this 
research. 
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IL  Pitch  Response  Criteria 


This  chapter  describes  the  aircraft  pitch  response  in  general  and  the  pitch  response  criteria 
examined  in  this  research.  First,  the  background  material  on  the  short  period  approximation 
is  given.  Next,  the  description  and  pertinent  information  on  the  Ralph  Smith  and  the  Time 
Domain  criteria  is  given.  Finally,  an  example  of  how  to  apply  the  two  criteria  to  an  aircraft 
transfer  function  is  presented. 

2. 1  Background 


In  aircraft  control  theory,  it  is  not  uncommon  to  have  feedback  control  systems  of  twen¬ 
tieth  order  or  more.  Writing  a  specification  for  such  a  large  system  can  be  cumbersome. 
Considerable  research  has  been  devoted  to  reducing  the  order  of  these  high-order  feedback 
control  systems  by  matching  frequency  responses  to  obtain  lower-order  equivalent  systems. 
Using  lower-order  equivalent  systems  allows  the  application  of  well-established  boundaries 
generated  by  classical  airplane  data  to  be  extended  to  many  high  order  systems  [1].  The 
pitch  angle  transfer  function,  £-,  for  a  linearized,  reduced-order  model  of  the  aircraft  is  given 


as 

9{s)  Kg  (Tg^  •  S  +  1)  {Tg^  •  S  -f  1)  6 

where 

^e{s)  +  2CpU}p  •  S  -h  (§2  +  2C,pWsp  ’  «  + 

e 

Pitch  Angle 

<5e  — 

Elevator  Deflection 

Kg  — 

Pitch  Transfer  Function  Gain 

T.,  — 

Low  Frequency  Pitch  Attitude 

Tfla - 

High  Frequency  Pitch  Attitude  Zero 

Tg - 

Aircraft  Time  Delay 

Cp 

Damping  Ratio  of  the  Phugoid  Mode 

LOp 

Undamped  Natural  Frequency  of  the  Phugoid  Mode 

Csp 

Damping  Ratio  of  the  Short  Period  Mode 

^sp 

Undamped  Natural  Frequency  of  the  Short  Period  Mode 

s  - 

Laplace  Variable 

Kq  can  be  further  defined  in  terms  of  the  non-dimensional  stability  derivatives  as 


Ke 


qiS‘^c[Cms^  {Clq,  +  Cp^)  -  CrrigCLs^ 
Iyy'^U\ 


(2.2) 
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where 


Qi  —  Dynamic  Pressure 

S  —  Wing  Reference  Area 

c  —  Mean  Aerodynamic  Cord 

Cma  —  Coefficient  of  Moment  due  to  Angle  of  Attack 

Cls^  —  Coefficient  of  Lift  due  to  Elevator  Deflection 

—  Coefficient  of  Moment  due  to  Elevator  Deflection 
Gl^  —  Coefficient  of  Lift  due  to  Angle  of  Attack 
Cz)i  —  Coefficient  of  Drag  at  Equilibrium 
lyy  —  Moment  of  Inertia 

m  —  Mass  of  Aircraft 

Ui  —  Equilibrium  Velocity 

The  detailed  derivation  of  Eq.(2.2)  is  shown  in  Appendix  A. 

In  cases  where  the  forward  velocity  response  is  small  at  the  natural  frequency  of  the 
short-period,  further  reduction  is  possible  [12].  The  low  frequency,  low  damping  pole,  or 
phugoid  mode,  can  be  separated  from  the  high  frequency,  high  damping  short  period  mode. 
This  is  done  by  setting  the  forward  speed,  u,  to  zero  in  the  equations  of  motion  and  neglecting 
the  forces  in  the  X  direction,  since  they  contribute  mostly  to  the  changes  in  forward  speed. 
In  such  cases. 


^(s)  ^  Ke{Te,-s  +  l)e--o-^ 

<5e(s)  5  (s2  ^  2CspUJsp  '  S  -f 

may  be  used  in  place  of  Eq.(2.1)  and  is  called  the  short-period  approximation. 
For  example,  the  following  pitch  transfer  function, 

0(s)  -1.31(s  +  0.016)(s  +  0.3)e-^^ 

6i{s)  ~  (s2  +  0.00466s  0.0053)(s2  +  0.806s  +  1.311) 


(2.3) 


(2.4) 


is  for  an  four  engine  jet  transport  flying  straight  and  level  at  Mach  number  0.62  [12].  The 
short-period  approximation  for  Eq.(2.4)  is 


9{s)  __  -1.39(s-t-0.306)e-i® 
“  s(s2-f  0.805s  4-  1.325) 


Fig. (2.1)  compares  Bode  plots  of  the  original  higher  order  system  to  the  short-period  ap¬ 
proximation.  For  this  example,  the  short-period  approximation  describes  the  aircraft  pitch 
response  fairly  well  for  the  frequency  range  given. 
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2.^  Ralph  Smith  Handling  Qualities  Criterion 

The  Ralph  Smith  criterion  combines  time  response  methods  of  determining  handling 
qualities  with  frequency  response  methods  [3].  Time  response  methods  are  those  that  use 
the  response  to  an  input  to  relate  handling  qualities  with  aircraft  parameters  such  as  rise 
time,  or  settling  time.  Frequency  response  criteria,  on  the  other  hand,  predict  the  handling 
qualities  by  relating  parameters  to  a  pilot  model.  The  Ralph  Smith  criterion  is  an  open-loop 
criterion,  as  shown  in  Fig.  (1.2),  but  was  derived  using  an  optimal  pilot  model  as  well  as 
flight  test  data  from  Neal  and  Smith  [13]. 

The  Ralph  Smith  criterion,  as  well  as  the  Time  Domain  criterion,  is  three-dimensional, 
as  seen  in  Table  2.1.  The  criterion  consists  of  three  parameters,  time  to  first  peak,  tq, 
average  slope,  S,  and  phase  lag,  In  order  to  determine  the  handling  qualities 

level,  one  must  go  through  all  three  parameters.  The  overall  value  is  determined  by  the 
worst  rating  of  the  three  parameters.  For  example,  if  the  three  ratings  are  1,2,  and  3,  then 
the  overall  handling  qualities  level  would  be  Level  3. 
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Table  2.1:  Ralph  Smith  Criterion 


Parameter 

Level  1 

Level  2 

Level  3 

tq  (secs) 

None 

None 

<  -2 

None 

None 

<  -123° 

<  -165° 

<  -180° 

When  a  pilot  is  determining  the  handling  qualities  of  an  aircraft,  there  are  many  ques¬ 
tions  that  must  be  answered.  Two  of  the  important  of  these  questions  are  as  follows: 

Is  the  response  too  abrupt  or  too  sluggish? 

Does  the  response  require  considerable  pilot  compensation?  [3] 


The  first  question  can  be  answered  by  time  history  analysis,  while  the  second  question  can 
be  answered  by  looking  at  the  frequency  response  of  the  aircraft. 

2.2.1  Time  to  First  Peak,  tq 

The  time  to  first  peak  parameter,  tq,  is  defined  as  the  time  to  first  peak  of  the  pitch  rate 
response,  q{t),  to  a  step  input  of  stick  force.  If  the  response  is  over-damped,  tq  is  defined 
as  the  time  to  90  percent  of  the  steady  state  value[3].  The  lower  bound  is  an  approximate 
representation  of  the  limit  on  human  time  delay,  litq  <  .2  seconds  the  pilot  tends  to  chase 
the  response.  The  typical  pilot  comment  would  be  that  the  aircraft  response  is  too  abrupt. 
With  a  time  to  first  peak  less  than  .2  seconds,  precision  maneuvering  will  be  difficult  without 
excessive  pilot  compensation  [3].  The  upper  bound  is  set  from  the  Neal-Smith  flight  test 
data.  Prom  pilot  comments  it  was  noted  that  a  system  with  tq  >  .9  tends  to  be  too  sluggish. 
This  results  from  excessive  lag  in  the  phase  angle  of  -^{jujc)  [3].  Fig.  (2.5)  shows  a  plot  of 
the  average  Cooper-Harper  rating  given  to  the  Neal-Smith  flight  test  data  versus  tq.  The 
numerical  data  is  given  in  Appendix  B.  By  looking  at  this  data,  it  can  be  seen  that  these 
bounds  are  adequate,  since  no  Level  1  ratings  occur  outside  the  bounds  which  mark  the 
Level  1  region. 
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Figure  2.2:  Average  Cooper-Harper  Rating  vs.  Time  to  First  Peak  of  Neal-Smith  Test  Data 
2.2.2  Slope  Parameter,  S 

The  slope  parameter,  5,  is  defined  as  the  average  slope  of  the  magnitude  plot  of  the 
transfer  function  ^{joj)  on  the  frequency  range  of  1  <  a;  <  6  radians/second.  This  slope  is 
representative  of  the  sensitivity  of  the  response  to  pilot  technique  [3] .  The  parameter  takes 
into  account  the  variability  of  pilots  by  requiring  the  slope  to  be  small,  thus  making  the 
aircraft  resistant  to  different  pilot  techniques  or  skill  level.  The  magnitude  of  the  slope  can 
be  determined  using  a  least  squares  best  fit  straight  line  on  the  frequency  range. 

The  boundaries  were  determined  by  dividing  the  Neal-Smith  flight  test  data  into  three 
groups, 

S  >  —2dB/oct  (2.6) 

—2  >S  >  —6dB/oct  (2.7) 

S  <  —6dB/oct  (2.8) 

Eq.(2.6)  yields  degraded  handling  qualities  and  is  too  abrupt,  similar  to  a  pure  gain  forcing 
function,  Kc-  Eq.(2.7)  produces  good  results  if  ljr{ju;c)>  -130°.  If  this  is  true,  the  aircraft 
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performance  is  governed  by  the  forcing  function  This  forcing  function  represents  a  simple 
stereotype  of  the  general  classification  of  aircraft-FCS  dynamics  which  was  determined  in 
McRuer’s  experiments  in  [15].  ^  turned  out  to  have  the  optimum  handling  qualities  of 
all  the  forcing  functions  tested.  Eq.(2.8)  acts  like  ^  which  exhibits  excessive  phase  lag  and 
is  stable  only  with  pilot  compensation  [3].  The  criterion  ignores  the  region  described  by 
Eq.(2.8),  because  it  is  usually  not  possible  to  get  a  Level  1  aircraft,  and  sets  the  boundaries 
as  S'  <  —2  dB/octave.  The  results  are  shown  in  Fig. (2.3)  using  the  flight  test  data.  No 
Level  1  ratings  were  given  to  aircraft  with  S  >  — 2  dB/octave. 


Figure  2.3:  Average  Cooper-Harper  Rating  vs.  Slope  of  Neal-Smith  Test  Data 
2.2.3  Phase  Lag  Parameter, 

The  phase  lag  quantifies  the  level  of  pilot  compensation  needed  to  perform  maneuvers. 
The  criteria  levels  were  determined  from  flight  test  data  and  pilot  comments.  Physically, 
phase  lag  is  the  amount  of  time  between  the  input  of  a  command  and  when  the  response  of 
the  aircraft  is  noticed  by  the  pilot.  In  order  to  calculate  l^{juJc),  the  criterion  frequency,  Wc, 
needs  to  be  determined.  This  criterion  frequency  is  approximately  the  crossover  frequency  of 
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the  pilot-aircraft  system  for  pitch  angle  tracking.  It  was  determined  by  using  the  crossover 
frequency  of  the  forcing  functions,  KcKJs,  and  Kc/s"^  from  McRuer’s  experiments  [15]. 
By  plotting  these  crossover  frequencies  against  the  forcing  functions’  slopes  in  dB/octave, 
the  criterion  frequency  can  be  defined.  Fig. (2.4)  shows  that  the  criterion  frequency  is  given 
by  the  equation  of  the  best  fit  straight  line  through  the  crossover  frequencies. 

Uc  =  .245  -b  6.0  (2.9) 

Once  the  criterion  frequency  is  calculated,  can  be  found.  This  is  done  by 

locating  the  phase  angle  at  on  the  Bode  phase  plot. 


Figure  2.4:  Specification  of  the  Criterion  Frequency 


2.3  Time  Domain  Handling  Qualities  Criteria 

The  Time  Domain  criterion  avoids  the  identification  of  dominant  roots  or  equivalent 
systems  models  by  working  directly  with  the  pitch  rate  transient  response.  Fig.  (2.5)  shows 
a  typical  pitch  rate  time  history.  The  following  measurements  are  defined  [16,  1]: 
a.  A  horizontal  line  defining  the  steady-state  pitch  rate,  Qss- 
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b.  A  sloping  straight  line  tangent  to  the  pitch  rate  time  history  at  the  point  of 
maximum  slope.  It  is  extended  to  intersect  both  the  steady  state  line  and  the 
time  axis. 


c.  Time  ti  measured  from  the  instant  the  step  input  is  applied  to  the  time  of 
intersection  of  the  maximum-slope  line  with  the  time  axis. 

d.  Time  t2  measured  from  the  instant  the  step  input  is  applied  to  the  time  cor¬ 
responding  to  the  intersection  of  the  maximum-slope  line  with  the  steady-state 
line. 


e.  Agi  =  maximum  pitch  rate  minus  the  steady  state  value 

f.  A^2  =  steady  state  minus  time  to  the  first  minimum. 


The  parameters  defined  above  should  meet  the  requirements  described  in  the  following 
subsections  when  a  step  input  elevator  deflection  or  a  step  stick  force  is  applied. 

2.3.1  Transient  Peak  Ratio,  Aq2/Aqi 

The  Time  Domain  criterion  are  stated  in  terms  of  the  transient  peak  ratio,  Aq2/Aqi  or 
TPR.  This  is  done  to  ensure  that  there  is  enough  damping  of  the  short  period  mode  of  the 
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pitch  response.  The  specific  values  are  based  on  the  interpretation  of  short-period  data  in 
[17]  and  [18]  and  are  shown  in  Table  2.2. 


Table  2.2:  Transient  Peak  Ratio,  Aq2/Aqi,  Levels 


Level 

Maximum  A52/A51 

1 

<  .30 

2 

<  .60 

3 

<  .85 

2.3.2  Equivalent  Time  Delay,  ti 

The  time  ti  is  considered  the  equivalent  time  delay  and  can  be  uniquely  defined  graphi¬ 
cally,  as  seen  in  Fig.  (2.5)  [16].  In  order  to  calculate  ti,  the  x-intercept  is  determined  from  the 
equation  of  the  line  tangent  to  the  maximum  slope  point.  The  limits  shown  in  Table  2.3  [16] 
were  determined  by  looking  at  flight  test  data  from  [19]  and  [20]  for  the  terminal  flight 
phase. 


Table  2.3:  Equivalent  Time  Delay,  ti.  Levels 


Level 

Equivalent  Time  Delay 

1 

ti  <  .12sec 

2 

ti  <  .17sec 

3 

ti  <  .21sec 

2.3.3  Rise  Time  Parameter,  At 

The  rise  time  parameter.  At,  is  defined  as 


At  =  t2  —ti. 


(2.10) 


Table  2.4  shows  the  requirements  that  must  be  met  for  the  rise  time  parameter  for 
terminal  flight  phases.  No  Level  3  is  defined.  These  limits  are  derived  directly  from  the 
limits  on  Un^/n/a  [16]. 


_  Qinitial  Qss/ At  _  Q 

~  Qss'^  ~  VxAt 


(2.11) 
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The  limits  on  /njoL  are  defined  as  a  function  of  the  different  flight  phases  [1].  Using  these 
constant  limits,  it  is  easy  to  show  that  the  boundaries  for  At  are  only  a  constant  divided 
by  the  true  airspeed.  For  example,  for  an  aircraft  in  landing  phase  the  Level  1  limits  are 

U)  ^ 

.16  <-^<3.6  (2.12) 

n/a 

Substituting  in  Eq.(2.11)  and  rearranging  for  At  gives 


<  At  < 


3.6Ur  “  “  .IWt 

Table  2.4  [16]  shows  the  limits  that  were  derived  from  the  data,  with  Vt  in  feet/second. 


(2.13) 


Table  2.4:  Rise  Time  Parameter,  At  Levels 


Level 

Minimum  At 

Maximum  At 

1 

9/Vt 

200/Ur 

2 

3.2jVT 

645/Ur 

2.4  Example  Point 

This  section  will  consider  one  point  and  show  how  to  determine  the  flying  qualities  level 
for  both  the  Ralph  Smith  criterion  and  the  Time  Domain  criterion.  The  sample  transfer 
function  is 

e  _  (2.8553s  +  2.445)e-^® 

Te  ~  s(s2  +  9.5178s  +  46.2182) 

Fig.(2.6)  and  Fig.(2.7)  show  the  pitch  rate  response  and  Bode  plots,  respectively  needed 
for  this  example.  First,  the  pitch  rate  response  is  generated  using  a  unit  step  input  and  time 
delay  of  .1  seconds.  Next,  a  line  representing  the  steady  state  value  is  drawn  on  the  pitch 
rate  history  plot.  Then  a  line  tangent  to  the  point  of  maximum  slope  is  plotted.  The  time, 
ti  is  calculated  by  measuring  the  difference  between  the  time  the  input  is  applied  and  the 
time  the  maximum  slope  line  crosses  the  x  axis.  In  this  example,  U  =  .1  seconds.  Recalling 
Table  2.3,  this  value  has  a  rating  of  1.  The  next  step  is  to  determine  the  time  the  maximum 
slope  line  crosses  the  steady  state  line.  Subtracting  U  from  this  value  gives  At  =  .0185. 
From  Table  2.4  and  Vt  =  65.23  meters/second  (214  feet/second),  this  corresponds  to  a 
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rating  of  2.  Now  determine  the  transient  peak  ratio  by  computing  Aqi,  the  maximum 
pitch  rate  minus  the  steady  state  value,  and  A52,  the  steady  state  value  minus  the  first 
minimum  value.  Dividing  A52  by  A^i  gives  the  transient  peak  ratio,  TPR  =  .0406,  and 
from  Table  2.2  this  has  a  rating  of  1.  The  overall  Time  Domain  handling  qualities  rating 
would  be  a  Level  2,  since  At  has  the  worst  rating  of  the  three  parameters.  Using  the  same 
pitch  rate  response,  determine  tq,  which  is  the  time  the  maximum  pitch  rate  occurs.  For  this 
example,  tq  =  .2843  seconds.  From  Table  2.1,  this  has  a  rating  of  1.  Looking  at  Fig.(2.7) 
determine  the  average  slope  of  the  magnitude  plot  on  the  frequency  range  of  1  <  a;  <  6 
radians/second.  This  is  determined  by  doing  a  least  squares  fit  of  the  magnitude  over  the 
frequency  range.  For  this  example,  S  =  —1.4  dB/octave,  which  corresponds  to  a  rating  of  2. 
Now  using  Eq.(2.9)  and  the  slope  parameter  from  above,  determine  the  criterion  frequency. 
For  this  example  Uc  =  5.66  radians/seconds.  Locate  this  frequency  on  the  phase  plot  and 
to  get  the  corresponding  phase  lag  angle,  ijrijuc)  =  -116.35°.  Once  again  looking  a 
Table  2.1,  this  corresponds  to  a  rating  of  1.  The  overall  Ralph  Smith  rating  would  be  a 
Level  2,  because  the  worst  parameter  rating  is  from  S.  Table  2.5  summarizes  the  numerical 
value  of  the  parameters  and  corresponding  ratings,  as  well  as  the  overall  handling  qualities 
level  according  to  each  of  the  criteria. 


Table  2.5:  Data  For  Example  Point 


Time  Domain  Criterion 

Parameter 

Value 

Rating 

ti 

.1001  sec 

1 

At 

.0185 

2 

Aq2/Aqi 

.0460 

1 

Overall  Level 

2 

Ralph  Smith  Criterion 

Parameter 

Value 

Rating 

tq 

.2843  sec 

1 

S 

-1.4  dB/oct 

2 

-116.35° 

1 

Overall  Level 

2 
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lagnitude,  dB 


III.  Results 


This  chapter  outlines  the  approach  that  was  used  to  map  the  criteria  into  each  other  and 
determine  the  areas  of  agreement  and  disagreement. 

3.1  Basic  Approach 

In  order  to  compare  the  Ralph  Smith  and  Time  Domain  criteria  it  is  necessary  to  look 
at  both  criteria  in  a  common  arena.  Since  both  criterion  require  a  pitch  transfer  function, 
the  first  step  is  to  calculate  the  short  period  approximation  of  the  f-  transfer  function. 
Recalling  Eq.(2.3),  it  can  be  seen  that  there  are  five  variables  that  need  to  be  chosen,  LOgp, 
(sp,  Tgj,  Kg,  and  rg.  Of  these  parameters,  three  are  set  by  the  aircraft  configuration.  This 
leaves  only  cogp  and  (sp  to  be  determined.  By  varying  these  two  parameters  a  grid  can  be 
set  up  and  each  criterion  examined.  Before  going  into  the  specific  approach,  there  are  two 
concerns  that  must  be  addressed.  First,  since  both  criteria  specify  boundaries  for  Level  3, 
it  is  likely  that  some  aircraft  configurations  will  fall  outside  of  this  level.  Points  that  are 
labeled  as  loss  of  control,  LOG,  are  those  which  do  not  fit  the  criteria  bounds. 

Another  likely  problem  that  can  occur  is  in  an  over-damped  system.  This  has  two 
effects  on  the  handling  qualities  level.  First,  an  over-damped  system  causes  a  change  in 
the  definition  of  tq,  from  the  Ralph  Smith  criterion.  This  definition  change  can  alter  the 
predicted  handling  qualities  level,  by  creating  a  discontinuity  in  the  Level  1  region.  The 
place  where  this  definition  change  occurs  is  labelled,  in  this  research,  as  the  jump  line. 
The  other  effect  that  an  over-damped  system  has  in  the  Time  Domain  criterion,  with  the 
transient  peak  ratio  parameter.  Recalling  Section  2.3.1,  TPR  is  used  to  ensure  that  there 
is  enough  damping  in  the  pitch  rate  response.  The  parameter  does  not  address  what  to  do 
when  the  system  is  sluggish,  or  over-damped.  Therefore,  in  this  research  TPR  is  neglected 
when  evaluating  an  over-damped  system  and  the  flying  qualities  level  is  determined  using 
the  other  two  parameters. 
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3.2  Test  Cases 

Table  3.1  shows  the  24  different  test  cases  which  were  examined  in  this  research.  Every 
four  cases  represent  a  different  aircraft  transfer  function  and  flight  condition.  Each  flight 
condition  is  examined  at  four  different  time  delays,  Tg. 


Table  3.1:  Test  Cases 


Case 

Height 

Velocity 

Kg 

rg 

Remarks 

m  (ft) 

m/s  (ft/s) 

s 

1 

304.8  (1000) 

65.23  (214) 

1.17 

2.44 

0 

Learjet 

Transfer  Function  I 

2 

n 

3 

.15 

4 

.2 

5 

0 

65.23  (214) 

1.41 

3.08 

0 

Learjet 

Transfer  Function  II 

6 

n 

7 

.15 

.2 

9 

0 

51.82  (170) 

1.32 

3.9897 

0 

Learjet 

Transfer  Function  III 

10 

.1 

11 

.15 

12 

.2 

13 

1524  (5000) 

77.11  (253) 

2.12 

.237 

0 

F-16 

Transfer  Function  I 

14 

n 

15 

.15 

16 

17 

304.8  (1000) 

65.23  (214) 

1.69 

.1144 

0 

F-16 

Transfer  Function  II 

18 

B 

19 

.15 

20 

.2 

21 

30.48  (100) 

64.2  (212) 

1.77 

.1431 

0 

F-16 

22 

.1 

Transfer  Function  III 

23 

.15 

24 

.2 

3.3  Step  By  Step  Process 

This  section  will  go  through  the  steps  required  to  create  a  mapping  of  each  of  the  criteria. 
An  example  is  provided  to  illustrate  the  algorithm.  The  example  corresponds  to  Case  2  from 
Table  3.1.  Results  of  all  the  cases  examined  are  located  in  Appendix  C. 
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Step  1:  Fix  aircraft  and  flight  phase.  This  determines  the  aircraft  classiflcation  and 
category.  For  this  example,  the  Calspan  Learjet  24,  a  Class  II  aircraft,  will  be  chosen.  The 
aircraft  will  be  in  a  Category  C  flight  phase,  i.e.  landing. 

Step  2:  Set  the  altitude,  h,  and  the  velocity,  Vt-  This  is  needed  to  find  the  non- 
dimensional  stability  derivatives  which  correspond  to  the  flight  condition.  The  example 
has  the  aircraft  at  a  height  of  304.8  meters  (1000  feet)  and  a  velocity  of  65.23  m/sec  (214 
ft/sec). 

Step  3:  Determine  and  Kq  from  the  non-dimensional  stability  derivatives  and  Eq.(2.2). 
These  two  variables  only  effect  the  numerator  of  the  transfer  function.  Therefore, 

Nl  =2.U{1. 17  s  +  1)  (3.1) 

Step  4-'  Fix  the  time  delay,  rg.  For  the  example,  Tg  =  .1  seconds. 

Step  5:  Determine  the  Ralph  Smith  criterion  for  a  grid  of  transfer  functions  by  varying 
(jJap  and  Csp-  This  research  used  a  grid  of  10,000  points  to  determine  the  criterion  as  was 
previously  discussed  in  the  example  point  in  Chapter  II.  Fig.(3.1)  shows  the  result  of  this 
step  for  the  example. 

Step  6:  Determine  the  Time  Domain  criterion  in  a  similar  manner,  as  seen  in  Fig.  (3.2). 

Step  7:  Cross  plot  each  level  to  show  the  areas  of  agreement  and  conflict.  Fig. (3.3)  shows 
the  comparison  of  the  Level  1  regions.  The  Level  2  regions  are  compared  in  Fig. (3.4).  The 
Level  3  comparison  is  shown  in  Fig.  (3.5),  while  Fig. (3.6)  compares  the  two  loss  of  control 
regions  for  this  example. 

The  computer  code  that  resulted  from  the  algorithm  was  written  for  MATLAB'^^  and  is 
shown  in  Appendix  D.  It  examines  the  two  criteria  for  a  grid  of  10,000  points. 
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Figure  3.1:  Ralph  Smith  Criterion  for  Learjet  Example 


Level  1 


Level  2 


Level  3 


Loss  of  Control 


Figure  3.2:  Time  Domain  Criterion  for  Learjet  Example 


IV.  Analysis 


One  objective  of  this  research  was  to  examine  trends  resulting  from  changing  parameters. 
This  chapter  will  examine  the  results  from  that  study. 

4-1  Ralph  Smith  Boundaries 

Fig.  (4.1)  shows  where  each  of  the  definitions  of  the  Ralph  Smith  criterion  take  effect.  It 
can  be  seen  that  only  Level  1  has  multiple  boundaries.  This  is  because  it  is  the  only  level 
with  definitions  in  all  three  parts  of  the  criterion. 


Figure  4.1:  Ralph  Smith  Boundaries 


The  jump  line,  as  seen  in  Fig.  (4.1),  is  the  boundary  which  separates  the  two  definitions 
of  the  time  to  first  peak,  tq.  The  value  of  tq  never  goes  to  any  particular  value,  but  jumps 
from  the  one  definition  to  the  next,  thus  changing  the  flying  qualities  level.  This  jump  line 
is  the  beginning  of  the  discontinuity  which  is  present  in  the  Level  1  region.  As  Csp  increases, 
the  pitch  rate  response  becomes  over-damped.  The  Ralph  Smith  criterion  states  that  for  an 
over-damped  system  tq  is  no  longer  defined  as  time  to  first  peak,  but  as  the  time  to  90%  of 
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the  steady  state  value.  Once  the  definition  changes,  the  pitch  rate  response  is  over-damped, 
and  tq,  using  the  new  definition,  becomes  smaller  than  the  minimum  bound.  Continuing  to 
increase  Csj?  increases  tq  so  it  is  within  the  Level  1  boundary,  temporarily.  Ultimately,  the 
pitch  rate  response  gets  too  slow  for  a  Level  1  rating  and  is  given  a  Level  2,  since  no  other 
limits  are  defined.  An  example  of  two  transfer  functions  on  either  side  of  the  jump  line  is 
shown  below. 


d{s)  (2.8533s -b  2.445)e--^" 
6e{s)  ~  s(s2-b  41.793s -1-36.238) 


and 

e{s)  (2.8533s -b  2.445)e-^* 

~  s(s2  -b  43.346s  -b  36.238) 

both  represent  aircraft  transfer  functions,  with  the  same  numerator  and  (Vgp-  The  only 
difference  is  in  the  short  period  damping  ratio:  =  3.4713  while  Csp2  =  3.6003.  Fig.  (4.2) 

and  Fig. (4.3)  show  the  pitch  rate  response  for  Eq.(4.1)  and  Eq.(4.2),  respectively.  By 
changing  Cspj  the  pitch  rate  response  goes  from  having  an  overshoot  and  a  rating  of  1,  to 
being  over-damped  with  a  rating  of  2. 


4.1.1  Effects  of  Changing  7$  on  the  Ralph  Smith  Criterion 

This  section  will  examine  the  effect  that  changing  the  time  delay,  r^,  has  on  the  Ralph  Smith 
criterion.  The  levels  will  be  broken  down  into  two  areas  of  discussion  The  first  one  describes 
the  transformation  of  the  Level  1  region  as  tq  is  increased,  while  the  other  discusses  the 
Level  2  and  3  regions. 

Fig.  (4.4)  shows  the  Ralph  Smith  criterion  Level  1  region  as  it  transforms  with  changing 
time  delay,  tq.  When  =  0  seconds  the  region  shows  that  the  previously  mentioned 
discontinuity  is  present,  but  not  completely  visible.  As  tq  is  increased  to  .1  seconds,  the 
discontinuity  is  larger  and  the  entire  level  1  region  is  moved  upward  and  to  the  right.  The 
lower  boundary,  set  by  lfr{juc)  rotates  slightly  counter-clockwise.  Increasing  rg  to  .15 
seconds  removes  the  discontinuity  completely.  The  reason  for  this  is  that  the  system  is  no 
longer  receiving  a  rating  of  1  in  the  phase  lag  parameter  when  the  definition  of  tq  changes, 
thus  there  is  no  longer  a  jump  line.  Once  the  jump  line  and  discontinuity  disappear,  only 
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Figure  4.3:  Pitch  Rate  Response  for  Eq.(4.2) 


4-3 


sp 


Figure  4.4:  Effect  of  Changing  To  on  Level  1 

three  boundaries  are  left  in  this  region.  When  tq  is  increased  to  .2  seconds,  the  region 
rotates  in  a  similar  fashion,  while  the  tq  =  .2  boundary  is  no  longer  present. 

Fig.  (4.5)  is  the  transformation  of  the  Level  2  and  Level  3  regions  as  a  result  of  changing 
Tg.  In  order  to  avoid  confusion  the  Level  1  region  is  omitted  in  since  the  boundaries  change 
in  the  same  manner  as  seen  in  Fig.  (4.4).  Both  the  lower  Level  2  boundary  and  the  Level  3 
boundaries  are  defined  by  one  parameter,  Since  a  time  delay  is  represented  on  a 

the  Bode  plots  as  — r^o;,  meaning  the  phase  angle  decreases  as  a  function  of  the  frequency,  u. 
This  phase  angle  is  added  to  the  angle  of  the  transfer  function  without  the  time  delay,  to  get 
the  total  phase  angle.  This  means  that  as  Tg  increa.ses,  the  phase  lag  parameter  decreases. 
The  more  lag  in  the  system,  the  worse  the  handling  qualities  become,  thus  moving  the 
boundaries  of  the  Level  2  and  Level  3  regions  for  the  Ralph  Smith  criterion. 

4.1.2  Effect  of  Changing  on  the  Ralph  Smith  Criterion 

For  Tg  —  .1  seconds,  as  Te2  increases,  the  lower  boundary  of  the  Level  3  region  changes, 
as  shown  in  Fig. (4.6).  The  bulge  in  the  region  moves  upward  and  increases  in  size  as  Tg^ 
increases. 
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Figure  4.6:  Effect  of  Changing  on  the  Ralph  Smith  Criterion 


4-2  Time  Domain  Boundaries 

Fig.  (4.7)  shows  where  each  of  the  definitions  of  the  Time  Domain  criterion  take  effect. 


^sp 

Figure  4.7:  Time  Domain  Boundaries 


4-2.1  Effects  of  Changing  tq  on  the  Time  Domain  Criterion 

The  effective  time  delay,  ti,  corresponds  to  the  time  delay,  tq,  by  the  following  relation¬ 
ship. 

Te  ^  ti  (4.3) 

This  corresponds  to  the  curved  boundaries  on  the  lower  left  of  both  level  2  and  level  3.  The 
effect  of  changing  the  time  delay  can  be  seen  in  Fig.(4.8).  Since  tq  emulates  ti,  when  the 
time  delay  is  increased  beyond  the  equivalent  time  delay  boundaries  defined  in  Table  2.3, 
the  rating  of  ti  changes.  This  in  turn  completely  eliminates  the  lower  handling  qualities 
level.  When  rg  <  .12  seconds  all  the  regions  are  present.  If  .12  <  rg  <  .17  seconds,  the 
entire  Level  1  region  has  disappeared,  as  shown  in  Fig.(4.8-B),  leaving  Level  2  ratings  in  its 
place.  At  Tg  >  .17  seconds,  the  Time  Domain  criterion  no  longer  gives  a  region  of  Level  2, 
Fig.(4.8-C).  Increasing  the  time  delay  greater  than  .21  seconds,  leaves  only  a  loss  of  control 
region,  since  the  Level  3  region  disappears. 
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Figure  4.8:  Effect  of  tq  on  Time  Domain  Criterion 
4.2.2  Effect  of  Changing  on  the  Time  Domain  Criterion 


If  To  is  held  constant  while  varying  the  Time  Domain  Level  1  and  Level  2  shift  the 
upper  and  lower  boundaries.  Prom  MIL-STD-1797A,  there  exists  a  relationship  between 
ujgp  and  At  [1]. 


CAP 


u> 


$p 


n/a  VrAt 

CAP,  the  Control  Anticipation  Parameter,  is  related  to  Ugp  by 


(4.4) 


CAP  «  (4.5) 

Vj, 

By  substituting  in  the  boundary  values  given  for  At,  the  boundaries  in  ujsp  can  be  approxi¬ 
mated. 

Rearranging  for  ujgp  in  terms  of  At, 


1 


Te,At 


(4.6) 
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Figure  4.9:  Effects  of  Changing  Tg^  on  the  Time  Domain  Criterion 

By  substituting  in  the  boundaries  for  At,  it  turns  out  that  the  upper  boundary  value  for 
At  is  the  lower  boundary  in  ujgp  and  vice  versa.  Using  this  relationship  it  is  possible  to  see 
that  when  Tg^  is  increased  the  boundaries  are  expanded,  as  shown  in  Fig.  (4.9). 

4.2.3  Relationship  between  TPR  and  ( 

From  MIL-STD-1797A,  the  relationship  between  transient  peak  ratio  and  (sp,  for  the 
classical  aircraft  response,  is  shown  in  Table  4.1  [1].  The  limits  on  TPR  correspond  to  the 
lower  bounds  on  (sp-  They  are  fixed  for  all  values  of  rg  and  Tg2- 


Table  4.1:  Corresponding  Damping  Ratios  for  Transient  Peak  Ratios 


Level 

Ac?2/Agi 

C 

1 

.30  j 

.36 

2 

.16 

3 

...  ^85 

.052 

3 


4.3.1  Level  1  Regions 

The  area  of  agreement  for  Level  1  is  dependent  upon  the  value  of  Tg.  The  only  agreement 
between  the  Ralph  Smith  criterion  and  the  Time  Domain  Criterion  occurs  when  =  0 
seconds.  This  region  of  agreement  can  change  with  Tg^  as  seen  in  Section  4.2.2.  As  Tg 
increases,  the  region  of  agreement  no  longer  exists.  There  are  two  reasons  explaining  this 
trend.  First,  as  described  in  Section  4.1.1,  the  Ralph  Smith  criterion  Level  1  region  shifts  as 
the  time  delay  increases.  The  effect  of  this  shifting  moves  the  two  regions  away  from  each 
other.  Next,  when  Tg  becomes  greater  than  .12  seconds,  there  is  no  longer  a  Time  Domain 
Level  1  and  therefore  no  agreement.  The  entire  transformation  of  the  Level  1  region  of 
agreement  is  shown  in  Fig. (4. 10)  for  increasing  Tg. 
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4-3.2  Level  2  Regions 

The  Level  2  region  of  agreement  can  consist  of  discontinuous  area,  depending  on  the 
value  of  tq.  When  tq  =  0,  the  regions  of  agreement  consist  of  four  distinct  areas.  This  is 
because  the  Level  1  regions  are  present  and  overlapping,  Fig.(4.11-A).  If  rg  =  .1  seconds, 
Fig.(4.11-B)  shows  that  there  are  two  regions  of  agreement.  This  again  is  due  to  the  presence 
of  the  Level  1  regions.  After  tq  becomes  greater  than  .12  seconds,  there  is  only  one  region  of 
agreement,  since  all  of  the  Time  Domain  Level  1  has  disappeared  seen  in  Fig.(4.11-C).  At 
Tg  >  .17  seconds  there  is  no  agreement  between  the  criteria  in  Level  2,  Fig.(4.11-D),  since 
the  Time  Domain  no  longer  has  a  Level  2  region. 

4-3.3  Level  3  Regions 

In  a  similar  manner  to  the  Level  2  region  of  agreement.  Level  3  has  a  discontinuous 
agreement  area.  Fig.  (4. 12)  shows  the  transformation  of  the  area  as  Tg  increases.  The  largest 
amount  of  agreement  between  the  criteria  exists  when  Tg  =  .2  seconds,  at  this  condition  all 
of  Ralph  Smith  Level  3  maps  into  Time  Domain  Level  3,  seen  in  Fig.(4.12-D). 


C.  "10  =.15  sec  D.  1^0  =.2  sec 


Figure  4.11:  Comparison  of  Level  2  Regions 
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Figure  4.12:  Comparison  of  Level  3  Regions 


4.3.4.  Loss  of  Control  Regions 


All  of  the  Time  Domain  criterion  that  is  rated  as  LOG  maps  directly  into  the  Ralph 
Smith  region  of  LOG.  Fig.  (4. 13)  is  representative  of  the  ihapping  for  the  loss  of  control 
region  at  all  rg. 
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V.  Conclusions  and  Recommendations 


5. 1  Conclusions 

All  the  objectives  for  this  research  were  accomplished.  First  a  computer  code  was  de¬ 
veloped  to  determine  the  handling  qualities  level  for  both  the  Ralph  Smith  criterion  and 
the  Time  Domain  criterion,  simultaneously.  The  code  allows  a  grid  of  transfer  functions  to 
be  analyzed  for  a  set  aircraft.  A  system  for  mapping  the  areas  of  agreement  and  disagree¬ 
ment  was  then  developed.  By  systematically  tagging  each  point  on  the  grid,  the  conflicting 
regions,  as  well  as  the  regions  of  agreement  could  be  identified. 

With  the  cases  examined  in  this  research,  a  general  comparison  between  the  Ralph  Smith 
criterion  and  the  Time  Domain  criterion  can  be  discussed.  The  best  comparison  of  the  two 
criteria  is  in  the  Loss  of  Control  region.  Here  all  points  defined  as  LOG  in  the  Time  Domain 
criterion  are  also  LOG  in  the  Ralph  Smith  criterion,  regardless  of  the  time  delay.  The  only 
other  time  that  all  points  of  one  criterion  mapped  directly  into  another  is  the  Level  3  region 
of  the  Ralph  Smith  criterion  when  the  time  delay  is  equal  to  .2  seconds.  This  is  because 
the  Time  Domain  criterion  only  consists  of  a  Level  3  region  and  a  LOG  region.  Looking  at 
the  rest  of  the  cases  show  that  there  is  little  agreement  between  the  two  criteria,  especially 
in  the  Level  1  region,  since  the  Time  Domain  criterion  has  no  Level  1  region  after  re  >  .12 
seconds.  Therefore,  the  analysis  of  the  results  show  that  the  Ralph  Smith  criterion  and  the 
Time  Domain  criterion  have  very  little  in  common. 

One  way  to  increase  the  area  of  agreement  between  these  two  criteria  is  to  modify  the 
Ralph  Smith  criterion  in  order  to  loosen  the  boundaries  on  the  phase  lag  parameter.  By 
doing  this,  the  Level  1  region  can  be  increased  to  include  more  area.  A  modification  that 
can  be  made  to  the  Time  Domain  criterion  is  to  change  the  rise  time  parameter  bounds. 
By  decreasing  the  minimum  boundaries  and  increasing  the  maximum  boundaries,  the  Level 
1  and  2  regions  can  be  increased  to  allow  more  comparability  between  the  aircraft. 
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5.2  Recommendations  for  Future  Research 

The  results  produced  by  this  research  lay  another  brick  in  the  foundation  started  by 
Kish  and  Kivioja,  though  this  is  not  the  end.  Further  research  can  be  done  on  the  different 
flight  phase  categories.  Selected  points  from  any  region  of  conflict  could  be  flown  in  a  flight 
test  program  to  determine  the  actual  handling  qualities.  Flight. testing  some  points  can 
determine  the  true  flying  qualities  level,  since  the  aircraft  is  a  nonlinear  system  with  a 
pilot  and  the  results  of  this  research  is  for  a  linear  approximation.  The  results  could  then 
be  analyzed  to  determine  whether  the  Ralph  Smith  or  Time  Domain  criteria  need  to  be 
refined,  combined,  or  even  completely  eliminated. 
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Appendix  A .  Derivation  of  Kg 


This  section  will  show  the  derivation  of  used  in  the  short  period  approximation  of  the 
pitch  transfer  function.  The  derivation  uses  the  dimensional  stability  derivatives,  which  are 
unique  for  aircraft  and  flight  conditions. 


9  _  +  Zs^Ma)s  +  [MctZs^  —  Mj^Zg) 

[/is(s^  —  +  Mg  +  Ma)s  +  —  Met)) 

Since  Zs^Ma  <  U^Mg^,  Zs^Ma  can  be  ignored.  Therefore, 

9  _  UiMs^s  —  Mg^Za  +  MetZs^ 

Se  Uis{s^  —  +  Mg  +  Ma)s  +  ( —  Ma)) 

Now  deflne 


^sp 


MgZot 

Ui 


-Met 


+  Mg  +  Mg) 


1  A  MaZs^  —  Ms^Za 

Tg^  UiMg^ 

Substituting  in  the  above  deflnitions, 

e  c/iM..(5  +  ^) 

Sq  Uis(^s^  T  cOgp) 

Rearranging  and  cancelling  Ui  gives 


9_  _  Mg^ _ [Tg^s  T  1) _ 

Tg^  s(s^  +  2CspW 

sp^  +  ^Ip) 

Therefore,  recalling  the  form  of  Equation  2.3, 


Or,  substituting  back  for  l/Tg^, 


Kg  = 


(A.1) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 


A-1 


K,= 


MgZs^  —  Ms^Z, 
Ui 


(A.9) 


The  dimensional  stability  derivatives  are  defined  in  terms  of  the  non-dimensional  stability 
derivatives  as  follows 


—  j 

^yy 

(A.IO) 

Mt,  = 

^yy 

(A.11) 

-q,S{Ci,  +  Cr,,) 

(A.12) 

m 

.  _-QiSC,s. 
m 

(A.13) 

Substituting  back  into  Equation  A.9 

QiScCms,  +  Cd,)  q.ScCm. 


K,= 


^yy 


m 


lyy  m 


(A.14) 


Or  rearranging 


Kf,  = 


+  Cp.)  -  C^,ClJ 

Iyy'^U\ 


(A.15) 
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Appendix  B.  Neal-Smith  Test  Data 


Table  B.l  is  data  from  Reference  [3,  13]  with  the  parameters  the  Ralph  Smith  criteria 


Table  B.l:  Neal-Smith  Test  Data 


Case 

tg 

S 

OJc 

Average 

Cooper-Harper 

Rating 

lA 

.56 

-9 

-136 

5 

IB 

.58 

-8 

4.1 

-132 

3.25 

1C 

.65 

-8.5 

4 

-148 

4.5 

ID 

1.0 

-11 

3.2 

-148 

4.125 

IE 

1.35 

-13 

-162 

6 

IF 

1.85 

-15 

2.2 

-180 

8 

IG 

1.5 

-15 

2.2 

-204 

8.5 

2A 

.26 

-1 

6 

-97 

4.25 

2B 

.35 

-2 

6 

-130 

5.125 

2C 

.24 

-3 

5.5 

-89 

3 

2D 

.3 

-4 

5.2 

-107 

2.66 

■MM 

.45 

-4 

5.2 

131 

4 

1M 

.6 

-6 

IQQI 

-143 

3 

2G 

.6 

-6 

Biai 

-168 

7 

2H 

.7 

-8 

■a 

-157 

5.5 

21 

.85 

-8 

IIB 

-179 

8 

2J 

3 

-9 

3.8 

-171 

6 

3A 

.15 

-1 

6 

-63 

4.25 

3B 

.25 

M 

6 

-90 

4.5 

3C 

.3 

IM 

5.5 

-107 

3.5 

3D 

.35 

-5 

4.9 

-123 

4 

3E 

3.5 

-6 

-137 

4 

4A 

2 

6 

-135 

5.25 

4B 

.45 

2 

6 

7 

4C 

.57 

-2 

5.7 

8 

4D 

.65 

-5 

4.9 

-168 

8 

4E 

.7 

-6 

4.60 

-173 

7 

5A 

.3 

5.5 

6 

-144 

6 

uses. 
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Appendix  C.  Test  Cases  and  Results 


This  appendix  includes  all  the  cases  tested  in  this  research  study.  Table  C.l  shows  all  of 
the  cases  studied. 


Table  C.l:  Test  Cases  Examined 


Case 

Height 

Velocity 

Toi 

Ke 

n 

Remarks 

m  (ft) 

m/s  (ft/s) 

s 

1 

304.8  (1000) 

65.23  (214) 

1.17 

2.44 

0 

Learjet 

Transfer  Function  I 

2 

.1 

3 

.15 

4 

.2 

5 

0 

65.23  (214) 

1.41 

3.08 

0 

Learjet 

Transfer  Function  II 

6 

B 

7 

.15 

8 

.2 

9 

0 

51.82  (170) 

1.32 

3.9897 

0 

Learjet 

Transfer  Function  III 

10 

.1 

11 

.15 

12 

.2 

13 

1524  (5000) 

77.11  (253) 

2.12 

.237 

0 

F-16 

Transfer  Function  I 

14 

.1 

15 

.15 

16 

.2 

17 

304.8  (1000) 

65.23  (214) 

1.69 

.1144 

0 

F-16 

Transfer  Function  II 

18 

B 

19 

.15 

20 

.2 

21 

30.48  (100) 

64.2  (212) 

1.77 

.1431 

0 

F-16 

Transfer  Function  III 

22 

B 

23 

.15 

24 

Each  case  represents  a  unique  flight  condition.  The  data  given  in  Table  C.l  contains 
all  the  pertinent  information  required  to  create  the  numerator  of  the  transfer  function.  All 
cases  are  presented  with  six  figures,  in  a  similar  manner  to  the  example  of  Chapter  III. 
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C.  1  Case  1 


Lear  Jet  Transfer  Function  I 
e{s)  _  2.44(1.17s4-l)e°^ 

S  (§2  +  2CspCOspS  +  UJ^p) 


C-2 


10 


CL, 

c/a 
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10  10  10 
C 

sp 

Figure  C.3:  Level  1  Regions  (Case  1  ) 
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Figure  C.4:  Level  2  Regions  (Case  1  ) 
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Figure  C.5:  Level  3  Regions  (Case  1  ) 
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Figure  C.6:  Loss  of  Control  Regions  (Case 


C.2  Case  2 


e{s)  _  2.44(1.17s  +  l)e--^^ 

S  +  2(sp<^spS  +  U^sp) 
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Figure  C.7:  Ralph  Smith  Criteria  (Case  2  ) 
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Figure  C.8;  Time  Domain  Criteria  (Case 
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Figure  C.ll:  Level  3  Regions  (Case  2  ) 
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Figure  C.12:  Loss  of  Control  Regions  (Case 


C.  3  Case  3 


e{s)  _  2.44(1.17s  +  l)e-^^^ 

<5(s)  S  (§2  +  2Cp(^spS  +  co^p) 
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Figure  C.13:  Ralph  Smith  Criteria  (Case  3  ) 
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Figure  C.14:  Time  Domain  Criteria  (Case 
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Figure  C.15;  Level  1  Regions  (Case  3  ) 
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Figure  C.16:  Level  2  Regions  (Case 
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Figure  C.17:  Level  3  Regions  (Case  3  ) 
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Figure  C.18:  Loss  of  Control  Regions  (Case  3  ) 
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C.4  Case  4 


e{s)  _  2.44(1.17s  +  l)e-^^ 
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Figure  C.19:  Ralph  Smith  Criteria  (Case  4  ) 
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Figure  C.20:  Time  Domain  Criteria  (Case  4  ) 
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Figure  C.21:  Level  1  Regions  (Case  4  ) 
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Figure  C.22:  Level  2  Regions  (Case  4  ) 
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Figure  C.23;  Level  3  Regions  (Case  4  ) 
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Figure  C.24:  Loss  of  Control  Regions  (Case 
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C.5  Case  5 


Lear  Jet  Transfer  Function  II 

e{s)  _  3.08(1.415  + l)e°^ 

^(s)  S  (s2  +  2CspCOspS  + 
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Figure  C.25:  Ralph  Smith  Criteria  (Case  5  ) 
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Figure  C.26:  Time  Domain  Criteria  Case  5 
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Figure  C.27:  Level  1  Regions  (Case  5  ) 
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Figure  C.28:  Level  2  Regions  (Case 
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Figure  C.29:  Level  3  Regions  (Case  5  ) 


Figure  C.30:  Loss  of  Control  Regions  (Case  5  ) 
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C.  6  Case  6 


e{s)  _  3.08(1.41g  +  l)e-^^ 
<5(s)  s  (s2  +  2CspU}spS  +  a;2p) 
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Figure  C.31:  Ralph  Smith  Criteria  (Case  6  ) 
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Figure  C.32:  Time  Domain  Criteria  (Case 
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Figure  C.35:  Level  3  Regions  (Case  6  ) 
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Figure  C.36:  Loss  of  Control  Regions  (Case  6  ) 
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Figure  C.37:  Ralph  Smith  Criteria  (Case  7  ) 


Level  1 


Level  2 


Level  3 


Loss  of  Control 


Figure  C.38:  Time  Domain  Criteria  (Case 
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Figure  C.39:  Level  1  Regions  (Case  7  ) 


Ralph  Smith 
Level  2 


Time  Domain 


Level  2 


Regions  of 
Agreement 

Level  2 


Figure  C.40:  Level  2  Regions  (Case 
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Figure  C.41:  Level  3  Regions  (Case  7  ) 
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Figure  C.42:  Loss  of  Control  Regions  (Case 
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Figure  C.43:  Ralph  Smith  Criteria  (Case  8  ) 


Figure  C.44:  Time  Domain  Criteria  (Case  8  ) 
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Figure  C.47:  Level  3  Regions  (Case  8  ) 


Ralph  Smith 
Loss  of  Control 


Time  Domain 
Loss  of  Control 


Regions  of 
Agreement 

Loss  of  Control 


Figure  C.48:  Loss  of  Control  Regions  (Case 
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Figure  C.49;  Ralph  Smith  Criteria  (Case  9  ) 
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Figure  C.50:  Time  Domain  Criteria  (Case 
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Figure  C.51;  Level  1  Regions  (Case  9  ) 


Ralph  Smith 
Level  2 


Time  Domain 


Level  2 


Regions  of 
Agreement 

Level  2 


Figure  C.52:  Level  2  Regions  (Case 
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Figure  C.53:  Level  3  Regions  (Case  9  ) 
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Figure  C.54:  Loss  of  Control  Regions  (Case  9  ) 
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Figure  C.55:  Ralph  Smith  Criteria  (Case  10  ) 
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Figure  C.56:  Time  Domain  Criteria  (Case  10  ) 
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Figure  C.57:  Level  1  Regions  (Case  10  ) 
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Figure  C.58:  Level  2  Regions  (Case  10  ) 
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Figure  C.59:  Level  3  Regions  (Case  10  ) 
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Figure  C.60:  Loss  of  Control  Regions  (Case  10  ) 


C.ll  Case  1 1 


e{s)  _  3.9897(1.32s  + 

<5(s)  S  (s2  +  2CpOJspS  +  co^p) 


C-42 


Level  1 


Level  2 


Level  3 


Loss  of  Control 


Figure  C.61:  Ralph  Smith  Criteria  (Case  11  ) 
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Figure  C.62:  Time  Domain  Criteria  (Case  11  ) 
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Figure  C.63:  Level  1  Regions  (Case  11  ) 
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Figure  C.64:  Level  2  Regions  (Case  11  ) 


Ralph  Smith 
Level  3 


Time  Domain 


Level  3 


Regions  of 
Agreement 

Level  3 


Figure  C.65:  Level  3  Regions  (Case  11  ) 
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Figure  C.66:  Loss  of  Control  Regions  (Case  11  ) 
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Figure  C.67:  Ralph  Smith  Criteria  (Case  12  ) 
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Figure  C.68:  Time  Domain  Criteria  (Case  12  ) 
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Figure  C.69:  Level  1  Regions  (Case  12  ) 
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Figure  C.70:  Level  2  Regions  (Case  12  ) 
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Figure  C.71:  Level  3  Regions  (Case  12  ) 
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Figure  C.72:  Loss  of  Control  Regions  (Case  12  ) 
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Figure  C.73:  Ralph  Smith  Criteria  (Case  13  ) 
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Figure  C.74:  Time  Domain  Criteria  (Case  13  ) 
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Figure  C.75:  Level  1  Regions  (Case  13  ) 
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Figure  C.76:  Level  2  Regions  (Case  13  ) 
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Figure  C.77:  Level  3  Regions  (Case  13  ) 
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Figure  C.78:  Loss  of  Control  Region  (Case  13  ) 
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Figure  C.79:  Ralph  Smith  Criteria  (Case  14  ) 
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Figure  C.80:  Time  Domain  Criteria  (Case  14  ) 
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Figure  C.83:  Level  3  Regions  (Case  14  ) 
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Figure  C.84:  Loss  of  Control  Regions  (Case  14  ) 
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Figure  C.85:  Ralph  Smith  Criteria  (Case  15  ) 
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Figure  C.86:  Time  Domain  Criteria  (Case  15  ) 
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Figure  C.87:  Level  1  Regions  (Case  15  ) 
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Figure  C.88;  Level  2  Regions  (Case  15  ) 
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Figure  C.89:  Level  3  Regions  (Case  15  ) 
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Figure  C.90:  Loss  of  Control  Regions  (Case  15  ) 
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Figure  C.91:  Ralph  Smith  Criteria  (Case  16  ) 
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Figure  C.92:  Time  Domain  Criteria  (Case  16  ) 
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Figure  C.93:  Level  1  Regions  (Case  16  ) 


Ralph  Smith 
Level  2 


Time  Domain 


Level  2 


Regions  of 
Agreement 

Level  2 


Figure  C.94:  Level  2  Regions  (Case  16  ) 
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Figure  C.95:  Level  3  Regions  (Case  16  ) 
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Figure  C.96:  Loss  of  Control  Regions  (Case  16  ) 
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Figure  C.97:  Ralph  Smith  Criteria  (Case  17  ) 
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Figure  C.98:  Time  Domain  Criteria  (Case  17  ) 
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Figure  C.99:  Level  1  Regions  (Case  17  ) 
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Figure  C.lOO:  Level  2  Regions  (Case  17  ) 
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Figure  C.lOl:  Level 


3  Regions  (Case  17  ) 
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C.102:  Loss  ofContro)Regi„„3  (Case  17) 
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Figure  C.103:  Ralph  Smith  Criteria  (Case  18  ) 


Level  1 


Level  2 


Level  3 


Loss  of  Control 


Figure  C.104:  Time  Domain  Criteria  (Case  18  ) 
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Figure  C.107:  Level  3  Regions  (Case  18  ) 
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Figure  C.108:  Loss  of  Control  Regions  (Case  18  ) 
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Figure  C.109:  Ralph  Smith  Criteria  (Case  19  ) 
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Figure  C.llO:  Time  Domain  Criteria  (Case  19  ) 
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Figure  C.lll:  Level  1  Regions  (Case  19  ) 
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Figure  C.112:  Level  2  Regions  (Case 
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Figure  C.113;  Level  3  Regions  (Case  19  ) 
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Figure  C.114:  Loss  of  Control  Regions  (Case  19  ) 
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Figure  C.115:  Ralph  Smith  Criteria  (Case  20  ) 
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Figure  C.116:  Time  Domain  Criteria  (Case  20  ) 
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Figure  C.117:  Level  1  Regions  (Case  20  ) 
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Figure  C.118:  Level  2  Regions  (Case  20  ) 
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Figure  C.119:  Level  3  Regions  (Case  20  ) 
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Figure  C.120:  Loss  of  Control  Regions  (Case  20  ) 
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Figure  C.121:  Ralph  Smith  Criteria  (Case  21  ) 


Level  1 


Level  2 


Level  3 


Loss  of  Control 


Figure  C.122:  Time  Domain  Criteria  (Case  21  ) 
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Figure  C.123:  Level  1  Regions  (Case  21  ) 
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Figure  C.124:  Level  2  Regions  (Case  21  ) 
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Figure  C.127:  Ralph  Smith  Criteria  (Case  22  ) 
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Figure  C.128:  Time  Domain  Criteria  (Case  22  ) 


Ralph  Smith 
Level  3 


Time  Domain 


Level  3 


Regions  of 
Agreement 

Level  3 


Figure  C.131:  Level  3  Regions  (Case  22  ) 
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Figure  C.132:  Loss  of  Control  Regions  (Case  22  ) 
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Figure  C.133:  Ralph  Smith  Criteria  (Case  23  ) 
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Figure  C.134:  Time  Domain  Criteria  (Case  23  ) 


Ralph  Smith 
Level  1 


Time  Domain 


Level  1 


Regions  of 
Agreement 
Level  1 


Figure  C.135:  Level  1  Regions  (Case  23  ) 
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Figure  C.136:  Level  2  Regions  (Case  23  ) 


Ralph  Smith 
Level  3 


Time  Domain 


Level  3 


Regions  of 
Agreement 

Level  3 


Figure  C.137;  Level  3  Regions  (Case  23  ) 
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Figure  C.138:  Loss  of  Control  Regions  (Case  23  ) 
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Figure  C.139:  Ralph  Smith  Criteria  (Case  24  ) 
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Figure  C.140:  Time  Domain  Criteria  (Case  24  ) 
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Figure  C.141:  Level  1  Regions  (Case  24  ) 
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Figure  C.142:  Level  2  Regions  (Case  24  ) 
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Figure  C.143:  Level  3  Regions  (Case  24  ) 
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Figure  C.144:  Loss  of  Control  Regions  (Case  24  ) 


Appendix  D.  Computer  Code 


The  following  computer  code  was  written  for  MATLAB^^.  It  consists  of  two  script 
files.  The  first  one  produced  the  data  for  the  mappings  while  the  second  plots  the  data  to 
produce  scatter  plots  similar  to  those  in  Appendix  C. 

generate. m 

Vt=214;  '/.velocity  in  ft/sec 
Ttheta2=1.69; 
ktheta=. 1144; 
tau=. 15; 

'/.Set  up  grid  for  examining  criterion 
zeta=logspace(-l , 1 , 100) ; 
wn=logspace (-1,2,100); 

wnTtheta2=wn . *Ttheta2 ; 

overall=zeros(length(wn)*length(zeta) ,1) ; 

kount=l ; 

for  i=l : length(zeta) ; 
for  j=l : length(wn) ; 


'/.create  transfer  function  theta/deltae 

num=ktheta* [Ttheta2  1] ; 

den=[l  2*zeta(i)*wn( j )  wn(j)~2  0] ; 

'/.find  transfer  function  q/deltae 
[nuraq , denq] =minreal ( [num  0] , den) ; 

D-1 


numq=[numq(2)  numq(3)] ; 


“/.create  state  space  realization 
[a,b,c,d]=tf2ss(numq,denq) ; 

’/,  the  following  lines  come  from  step.m 

“/.  The  next  two  constants  control  the  precision  of  the  plot 

“/,  cind  the  time  interval  of  the  plot. 

st=0.005;  \*/.  Set  settling  time  bound  =  0.5“/. 

precision=30;  \“/.  Show  approx  30  points  for  simple  graph 

“/.  Step  response  is  effectively  equal  to  placing  initial  conditions 

“/.  on  the  plant  as  follows : 

[n,ra]=size(b) ; 
if  abs(rcond(a))  >  eps 
xO  =  -a\(b*ones(m, 1)) ; 

“/.  Cater  for  pure  integrator  case 
else 

xO  =  ones (n, 1) ; 
end 

t=t imvec (a, b,c,x0,st .precision) ; 

“/. 

“/.number  of  points  wanted  for  precision 

pts=10000; 

tconst=max(t) ; 

tconst=tconst+tau; 

“/.creates  em  ending  point  that  is  displaced  by  the  time  delay 

“/.set  up  a  linear  time  vector  which  starts  at  zero  and  ends  at  tconst  with 

“/.specified  number  of  points 

t=linspace(0, tconst ,pts) ; 
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y,add  time  delay  into  step  input 
u=stepfun(t ,tau) ; 

y, simulate  the  response 
q=lsim(numq,denq,u,t) ; 
y.plot  (t ,  q) 

‘/ofind  the  steady  state  value  of  the  response 
dcgain=-c/a*b+d ; 


y.find  qdot/deltae  transfer  function  inorder  to  find  slope 
numl= [numq  0] ; 
denl=[denq]  ; 

qdot=lsim(numl ,denl ,u,t) ; 

'/.find  time  that  max  slope  occurrs,  t(k)  and  point  on  q  response 
[maxslope,k] =max(qdot) ; 
qmaxslope=tablel( [t ’  q],t(k)); 

m=maxslope;  %  slope  of  maxslope  line 

‘/.find  the  time  tl  where  maxslope  line  crosses  y=0 
%y=0,  x=tl; 

b=(qmaxslope-maxslope*t(k)) ; 
tl=-b/m; 


’/.find  level  for  equivalent  time  delay 
if  tl<=  .12 
delaylevel(kount)=l ; 
elseif  tl<=  .17 
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delaylevel (kount ) =2 ; 
elseif  tl<=.21 
delaylevel (kount ) =3 ; 
else 

delaylevel(kount)=10; 

end 

7, find  time  t2,  where  maxslope  line  crosses  qss 
7oy=dcgain  x=t2 
t2= (dcgain-b) /m ; 

7, rise  time  parameter,  terminal  flight  phases 
deltat=t2-tl; 

if  deltat>=  9/Vt  &  deltat<=  200/Vt 

risetimelevel (kount) =1 ; 

elseif  deltat>=  3.2/Vt  &  deltat<=  645/Vt 

risetimelevel (kount) =2 ; 

else 

risetimelevel (kount ) =3 ; 
end 

7ofind  ql,  max  pitch  rate 
[ql,kl]=max(q) ; 

7ofind  time  to  first  peak 

if  ql  ==  q(length(q))  7olast  q  (System  is  overdamped) 

q90=.90*dcgain;  7ofor  overdamped  system  tq  is  defined  as 

Z=f ind(q>q90)  ;  7o.  time  to  907.  of  final  value 

tq=t(Z(l)); 
else 
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tq=t(kl) ; 
end 

'/.find  time  to  first  peak  level 
if  tq  >=  .2  &  tq  <=  .9 
tfplevel(kount)=l ; 
else 

tfplevel(kount)=2 ;  "/.no  other  level  are  defined 
end 

'/.find  deltaql 
deltaql=ql-dcgain ; 

'/.find  q2,  first  min 
[q2,k2]=min(q(kl:pts)) ; 

y.find  deltaq2 
deltaq2=dcgain-q2 ; 

'/.allow  for  possibility  of  no  min 
if  deltaq2<=  0 
deltaq2  =  0; 
end 

’/.allow  for  overdamped  system 

if  deltaq2~=0  &  deltaql<=0 

deltaq2=0 ; 

deltaql=0; 

end 
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'/.Transient  Peak  ratio 
tran=deltaq2/deltaql ; 
if  tran<=  .3 
tprlevel(kount)=l ; 
elseif  trcLn<=.6 
tprlevel (kount ) =2 ; 
elseif  trcin<=.85 
tprlevel (kount ) =3 ; 
elseif  tran>.85 
tprlevel (kount )=  10; 
else 
tran=0 ; 

tprlevel (kount) =1 ; 
end 


'/.produce  the  bode  magnitude  and  phase  plots  for  the  transfer  function 
w=logspace(0,logl0(6') ,  10) ; 

[magn , phase] =bode (num , den , w) ; 
mag=20*logl0(magn) ;  '/.put  mag  in  dB 

'/.find  slope  on  interval  of  l<w<6 
A=[logl0(w)’  ones(10,l)]; 
slope=inv(A'*A)*(A’*mag) ; 

S=.3*slope(l) ;  '/.-20db/dec=-6db/oct 

'/.Determine  the  Level  based  on  slope 
if  S  <  -2 

slopelevel(kount)=l ; 
else 
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slopelevel (kount ) =2 ; 
end 

'/.Find  criterion  Frequency 
wc=.24*S+6; 


“/oCreate  larger  interval  for  bode  plot 
wl=logspace (-2 , 2 , 300) ; 

[magi, phasel]=bode(num, den, wl) ; 
mag2=20*logl0(magl) ; 


Phase=phasel-tau*wl  ^  *180/pi ;  "/Adds  in  Time  delay 

’/.Determine  the  phase  angle  at  the  criterion  frequency 
lag=tablel( [wl’  Phase] ,wc) ; 

'/.Find  Phase  Lag  level 

if  lag  >  -123 

laglevel (kount )=1 ; 

elseif  lag  <=  -123  &  lag  >  -165 

laglevel (kount ) =2 ; 

elseif  lag  <=  -165  &lag  >  -180 

laglevel (kount ) =3 ; 

else 

laglevel(kount)=10; 

end 


'/.  overall  flying  qualities  level 
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timedomain(kount)=max( [delaylevel(kount)  risetimelevel(kount)  tprlevel(kouiit)] ) ; 
RS(kount)=inax( [tfplevel(kount)  slopelevel(kount)  laglevel(kount)] ) ; 

‘/.create  data  matrices 
'/.Ralph  Smith  Data 

RSdataCkount, :)=[zeta(i)  wnTtheta2(j)  RS(kount)  tfplevel(kount)  ... 
slopelevel(kount)  laglevel(kount)] ; 

'/.Time  Domain  Data 

TDdataCkount , : )= [zeta(i)  wnTtheta2(j)  timedomain(kount)  delaylevel(kount)  ... 
risetimelevel(kount)  tprlevel(kount)] ; 

'/.Data  for  both  criteria 

rawdataCkount, :)=[zeta(i)  wnTtheta2(j)  tl  Vt*deltat  tran  timedomain(kount)  ... 
tq  S  lag  RS(kount)  wn(j)] ; 

kount=kount+l ; 

end 

end 

save  /tmp_mnt/home/dynamics/lcarlucc/thesis/cases/f 16_2.mat 

exit 
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dataplot.m 


y  0/0/ 0/0/ 0/0/0/ 

/o  /e  /o  /o  /o  /o  /o  /o 

y.Maps  Time  Domain 

0/0/ 0/0/ 0/0/ 0/0/ 

/o  /o  /•  /9  /«  /o  /o  /o 

count=l ; 
figure 

loglog(0,0, ’k. ’) 

xlabtexC [gtex( ’z' )  ,stex('spO]) ; 
ylabtexC  [gtex(  'w’ )  ,stex(  'spOl ) 
hold 

for  i=l : length(zeta) ; 

for  j=l :length(wn) ; 
if  timedomainC count)  ==1 
loglog(zeta(i) ,wn(j) , 'rx’ ) 
elseif  timedomain(count)  ==  2 
loglog(zeta(i) ,wn(j) , 'yx’ ) 
elseif  timedomainC count)  ==  3 
loglogCzetaCi) ,wn(j) , 'bx') 
end 

count=count+l ; 

end 

end 


0/0/ 0/0/ 0/0/ 0/0/ 

/o  /q  /o  /o  /«  /o  /o  /o 

'/.maps  Ralph  Smith 

0/0/ 0/0/ 0/0/ 0/0/ 

/e  /o  /o  /o  /q  /o  /o  /o 

count=l ; 
figure 

loglog(0,0, ’k. ’) 


xlabtexC [gtex( 'z’ ) ,stex(’sp’)]) ; 
ylabtex ( [gtex ( ’ w ’ ) , stex ( ' sp ’ ) ] ) 
hold 

for  i=l:length(zeta) ; 

for  j=l : length (wn) ; 
if  RS(count)  ==1 
loglog(zeta(i) ,wn(j) , ’rx’ ) 
elseif  RS (count)  ==  2 
loglog(zeta(i) ,wn(j) , ’yx’) 
elseif  RS (count)  ==  3 
loglog(zeta(i)  ,wn(j) ,  ’bxO 
end 

count=count+l ; 

end 

end 

VVVVVVV 

/o  /o  /o  h  /o  /o  /« 

'/.Maps  level  1  of  TD,  RS 

v'/vvvvy 

h  /o  /o  /o  /e  /o  /e 

count=l; 

figure 

loglog(0,0, ’k. 0 

xlabtex ( [gtex ( ’ z  O , stex ( ^  sp  O ] ) ; 
ylabtex( [gtex( ’w’ ) ,stex(’sp’)]) 
hold 

for  i=l:length(zeta) 
for  j=l : length (wn) ; 

if  timedomain( count)  ==1 
loglog(zeta(i) ,wn(j) , 'bo') 
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end 

if  RS (count)  ==1 
loglog(zeta(i) ,wn(j) , ’yx’) 
end 

count=count+l ; 

end 

end 


vvvvvvv 

/o  /o  /«  h  /o  /e  /o 

'/.Maps  level  2  of  RS,TD 

VVVVVVy 

count =1 ; 
figure 

loglog(0,0, ’k. ’) 

xlabtexC [gtex( ’z’ ) ,stex( ’sp’ )] ) ; 
ylabtex([gtex(’w’) ,stex(’sp’)]) 
hold 

for  i=l : length(zeta) 
for  j=l : length(wn) ; 

if  timedomainC count)  ==2 
loglog(zeta(i) , wn( j ) , ’bo ’ ) 
end 

if  RS (count) ==2 
loglog(zeta(i) ,wn(j) , ’yx’) 
end 

count=count+l; 

end 
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end 


I 


VVV'/VVV 

/o/d/o/o/o/o/o 

y,Maps  level  3  of  TD,  RS 

vvvvvvv 

count=l ; 

figure 

loglogCO.O.’k.’) 

xlabtexC [gtex( ’z’) ,stex( ’sp’ )] ) ; 
ylabtex ( [gtex ( ’ w  O , stex ( ’ sp ’ ) ] ) 
hold 

for  i=l ; length(zeta) 
for  j=l : length (wn) ; 
if  timedoraainC count)  ==3 
loglog(zeta(i) ,wn(j) , 'bo’) 
end 

if  RS( count )==3 
loglog(zeta(i) ,wn(j) , 'yx') 
end 

count=count+l ; 

end 

end 

VVVVVVV 

/o  /o  /o  /o  /o  /o  /o 

‘/.Maps  level  10  of  TD,  RS 
VVVVVVV 

/o/o/o/o/o/o/o 

count=l ; 
figure 

loglog(0,0, 'k. ’) 
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xlabtexC [gtex( ’z’ ) ,stex( 'sp' )] ) ; 
ylabtex ( [gtex ( ’ w ’ ) , stex ( ' sp ’ ) ] ) 
hold 

for  i=l : length (zeta) 
for  j=l : length(wn) ; 
if  timedomainC count)  ==10 
loglog(zeta(i)  ,wn(j) ,  ’boO 
end 

if  RS (count) ==10 
loglog(zeta(i) ,wn(j) ,’yx’) 
end 

count=count+l ; 

end 

end 
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